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Abstract 

The slope of the Isgur-Wise function for ground state mesons is evaluated for 
the heavy mass limit of quark models a la Bakamjian-Thomas, which has been pre- 
viously discussed by us in general terms. A full calculation in various spectroscopic 
models with relativistic kinetic energy gives a rather stable result p 2 ~ 1, much 
lower than previous estimates. Attention is paid to a careful comparison of this 
result with the ones of QCD fundamental methods (lattice QCD, QCD sum rules) 
and with experimental data. 
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In a series of papers [JTJ| H [0] , we have shown that in the heavy mass limit tuq — > oo 
with one light quark, quark models of form factors based on the old Bakamjian-Thomas 
(B-T) construction of relativistic states, and the free quark current, are presenting very 
interesting features : first of all covariance ; but also heavy quark symmetry for form 
factors, i.e. heavy to heavy Isgur-Wise scaling behavior and normalisation [|1] ; duality 
properties amounting in particular to the Bjorken sum rule |§. And finally, we are showing 
in a forthcoming paper [|J that they satisfy the new sum rules involving annihilation 
constants that we have derived in ||. 

Now, of course, such general properties, so interesting as they may be, cannot be 
the last word. One has to tackle with concrete, quantitative predictions. This implies 
making a specific choice of the mass operator or spectroscopic model, which has been left 
unspecified up to now. It is this spectroscopic model which will determine the internal 
wave functions at rest of the B-T formalism and lead to the quantitative predictions that 
we begin to formulate in this letter. 

We will use actually various possible spectroscopic models that have been presented 
in the literature ; they are based mainly on the analysis of the spectrum, but indeed the 
spectrum is still roughly compatible with many possible models. In this situation, the 
study of transition rates can afford precious additional information and help selecting the 
spectroscopic models. At this point, one must warn the reader against a possible confusion 
in the comparison of models. In fact, several authors have proposed at the same time 1) a 
specific spectroscopic model and 2) a specific formulation of form factors in terms of wave 
functions, which differ from the B-T formulation ; but the two aspects must be discussed 
separately. In all the following, it should be clear that we are discussing the results of 
the various spectroscopic models of these authors within one and the same formulation 
of form factors, which is the B-T formulation advocated by us, except for some parts 
where we make also a comparison with the results obtained within the old standard non 
relativistic formulation of form factors. 

We begin with the calculation of a quantity which is the object of both a theoretical 
and experimental very active interest, the slope of the Isgur-Wise function at zero recoil, 
usually parametrized by : 



since t^-(I) < 0. The complete calculation of B — > D,D*,D* * lis transitions will be 
presented in a forthcoming paper ||. The main conclusion of the present paper is that a 
proper relativistic treatment of quark models (including the general Bakamjian-Thomas 
formulation and a relativistic spectroscopic model) gives values of p 2 around : 

p 2 ^l . (2) 

This value is considerably lower than a previous calculation by Close and Wambach 
0, who found p 2 = 1.4. While postponing detailed discussions to the end of the letter, 
we want to emphasize from the beginning that our number eq. (0) represents a very signi- 
ficant progress with respect to this previous higher estimate. Indeed, QCD fundamental 
methods, i.e. sum rules and numerical lattice QCD, are yielding rather low values. Sum 
rules predictions 

range from 0.54 to 1. UKQCD lattice calculation [12[ 
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yields 0.9 as central value, admittedly with very large error bars : 

p 2 = 0.9+g(stat)+2(syst). (3) 

Our quark model prediction p 2 pa 1 is therefore much closer to the predictions of 
QCD fundamental methods, than the 1.4 of ref. pj. 

The progress is even more striking if one considers s quarks instead of u quarks, 
because we find a larger difference with J7J, while UKQCD has naturally smaller errors. 
In this case, Close and Wambach found p 2 = 1.64, while our result is : 

^~1.15 . (4) 
The central value quoted by UKQCD is : 

^ = 1.2+2 (stat) +2 (syst) (5) 

It must be kept in mind, as we will explain later, that the results of quark models are 
not directly comparable to the renormalization-group invariant p 2 , which is the quantity 
given by QCD fundamental methods f], and which we henceforth denote as p 2 en , gr , inv ,- 
Nevertheless, the progress is quite encouraging. 

Of course, one would like also to confront theory to experiment, since the pre- 
dictions of QCD are not so safe and not so accurate. However, as we will explain, the 
message of experiment itself is not so clear. Roughly speaking, the experiment does not 
determine the slope at zero recoil as one commonly assumes, but an average slope over 
the allowed range of q 2 , which is sizeably different. Therefore, we prefer to postpone the 
discussion of experimental data until the end of the letter. 

It seems to us, meanwhile, already significant to get agreement with QCD methods, 
or at least to get much better compatibility with them than previously obtained in quark 
models. But then, how do we obtain this qualitative improvement ? There are two 
aspects. 

The main aspect is that we use relativistic spectroscopic models instead of a non- 
relativistic one in 0, i.e. models with a square root kinetic energy ^p 2 + m 2 instead of 
m + p 2 /2m. One must stress that the necessity of a relativistic expression of form factors 
(relativistic treatment of overall hadron motion and relativistic current operator) is quite 
obvious in calculating p 2 (otherwise one would get a very small result < 0.3); such a 
relativistic expression is indeed used by Close and Wambach 0j; but, while adopting such 
an expression, one can still think of maintaining a non-relativistic spectroscopic model 
to calculate wave functions at rest, because such models seem to describe the hadron 
spectrum reasonably this is what has been done by Close and Wambach 0. What 
we find however is that wave functions at rest can be strongly different in a relativistic 
spectroscopic model, although the spectrum is similar, whence the lowering of p 2 . 

Another, more technical, but important aspect is that the very common approxi- 
mation of using a variational Gaussian to approximate the wave function 0, |14|, fT5| , 



3 The first results along the same method were obtained by Bernard, Shen and Soni 

4 We read the number quoted by Narison in Jll[ as " the slope of " physical" Isgur-Wise function" as 

being actually the renormalisation-group invariant p 2 

5 This distinction between the general treatment of form factors and the spectroscopic models, each of 

them being possibly either relativistic or non-relativistic, will be present throughout the paper 
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fails even at low recoil, e.g. in the calculation of p 2 , especially in the context of 
relativistic spectroscopic equations. Calculating exactly the wave function results in an 
appreciably lower estimate of p 2 . 

To summarize, the improvement we obtain for p 2 is based on logical ingredients : 
more exact calculation of the wave function, and introduction of relativistic wave equa- 
tions, which seems unavoidable since the light quarks in a heavy light system are strongly 
relativistic. 



General statements on p 2 in the B-T approach 

The general expression of p 2 in our approach is given in ref. 
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In the paper ref. H] , an alternative formula is given, which expresses p 2 as a sum of three 
positive contributions, each having a simple meaning in a small velocity expansion of the 
matrix element < f'|7°|f > (normalisation 5 3 (p — p')) 
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(7) 



In both cases the matrix elements denoted by (0 |. . .| 0) are taken on ground state spatial 
internal wave functions ; p is the light quark momentum operator, and f is the canonically 
conjugate operator ; i.e., in momentum space which is the most convenient one for this 
calculation : 







(8) 



In eq. ([?]), the 1/4 corresponds to the contribution of the free active spin 1/2 heavy quark 
current to < i/|7°|i; > and it would be the same for any hadron when expressed in terms 
of the slope of this matrix element, instead of invariant form factors. The two other terms 
represent effects specific to composite hadrons, and depending on the considered hadron 
state. The first one would be present even if quarks would have zero spin. We denote it 
as p 2 space . The second one comes from the Wigner rotations of the spectator quark spin. 
We denote it as p 2 wigner . 

In [|TJ, we have found a lower bound on p 2 , larger than 1/4, and specific to models 
based on the B-T approach : 
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As explained in 
large \p\ 3> m is p 



this lower bound reflects the fact that the expression for p 2 „„„„ for 



space 



space 



X f \p\f+r\p\ 
3 V 2 



; then, from the uncertainty principle \p\r > 1, one 
naively expects p 2 space to be > 1. On the other hand for \p\ -C m (non-relativistic limit), 
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For a particular class of internal wave functions, the lower bound on p 2 must be 
obviously still larger. For Gaussian wave functions Jl| : 



Choice of the spectroscopic model 

As our preferred relativistic model to calculate the ground state internal wave func- 
tion <po(p 2 ), we choose the one of Godfrey and Isgur (GI) |l7fl. Certainly, this model is 
rather complex and because of that, one may not agree on all the ingredients or assump- 
tions which enter it, or one may find difficult to evaluate their respective impact on the 
final result. Nevertheless, what makes the model outstanding is its covering of the whole 
spectroscopy, from light to heavy quarks, 

Contrarily to other models in the heavy-light sector, the GI model contains the 
full set of spin-dependent forces. For the present purpose of calculating p 2 , we need the 
wave functions in the limit ttlq — ► oo. In this limit, spin-spin and tensor forces, as well 
as certain parts of the spin-orbit forces ( those coming from transverse gluon exchange), 
disappear. However, the full set of spin-dependent forces is required to describe the whole 
spectrum of states and therefore to assess the validity of the model. The GI model includes 
relativistic features, among which the root square kinetic energy : 



and the above various spin- dependent forces, but also many other effects correcting the 
gluon exchange potential by momentum dependent factors and smearing. Truely, the 
latter effects spoil the flavor independence of the potential since they depend on quark 
masses, and in a somewhat adhoc way. The kinetic energy eq. ([11]) is shared by several 
other models and we will also present the corresponding p 2 to settle the discussion. 

Program of numerical calculation of wave functions 

To calculate the eigenf unctions of the Godfrey-Isgur Hamiltonian, which will be 
used as the internal wave functions ip in the later calculation of p 2 , we use the method 
of these authors, i.e. to diagonalize the matrix of the Hamiltonian, taken on a large 
harmonic oscillator basis jl7| 0. We have made calculations with basis dimension N as 
large as 25. As one needs only the ijiq = oo limit, one can simplify much the program by 
dropping the spin-dependent contributions vanishing in this limit. A surprising feature of 
this limit is that the usual order of L = 1 levels is reversed, i.e. the 2 + is the lowest. We 
have tested the program by making rnq = m;, and comparing with the energy levels given 
in the paper of Godfrey and Isgur. The agreement is reasonable in view of the above 
simplifications. In fact, an exact check of the program has been obtained in the following 
way. One has also formulated the complete program which reproduces exactly the masses 
of Godfrey and Isgur at finite quark masses. Then, by dropping the relevant terms, one 
recovers exactly the numbers obtained with the simplified program. Programs and checks 
have also been made for the other spectroscopic models used in the discussions below. 

6 Similar calculations have been done to evaluate the pion and other light meson form factors 




(10) 




(11) 
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Results for p 2 and discussion of the wave functions 



Then p 2 is calculated straightforwardly. We find : 

p 2 GI = 1.023 . (12) 

As announced, this is much lower than the 1.4 result obtained in |7]], using a relativistic 
treatment of hadron motion, but a non-relativistic spectroscopic model for wave functions 
at rest, namely the ISGW spectroscopic model |14| 7 . In fact, the result of is obtained 



with a Gaussian approximation. Doing a complete numerical treatment for the same 
model, we obtain 1.283, lower than 0, but still much higher than our result eq. flT2| ) 
obtained with the GI spectroscopic model. 

On the contrary, for spectroscopic equations having a square root kinetic energy like 
the GI model, but differing by their potential part, as proposed for example by Veseli and 
Dunietz(VD) [19j which take a linear plus Coulomb potential, or by Colangelo, Nardulli 
and Pietroni(CNP) [p0[] , which use a Richardson potential, we find results remarkably 
similar to eq. flT2| ) : resp. p 2 = 0.98 ; 0.97 0. From the comparison of all these results, it 
seems justified to conclude that the decisive reason why we obtain a low value p 2 ~ 1, is 
the use of a relativistic kinetic energy instead of a non-relativistic one in [0. On the 
other hand, the result does not seem very sensitive to the very short distance behavior 
of the potential, since the three relativistic spectroscopic models give a similar p 2 , while 
their Coulomb-like potential a s (r)/r has respectively a regular short distance behavior 
for GI Q an 1/r singular behavior for since a s = est, a milder singularity for |HJ since 
a s ~ l/£nr. Neither is the result sensitive to the corresponding short distance behavior 
of (p, which is respectively regular, power-like singular, or logarithmically singular when 
r — ► [2T|. One must note, however, that, not unexpectedly, the convergence of the 



numerical method is slower for singular potentials with relativistic kinetic energy, and 
therefore we have less significant digits in this case. As announced, the new, low result 
p 2 « 1 cannot be obtained if one uses the very common approximation by a Gaussian 
wave function : 

(3 2 r 2 

(p(r) oc exp( — ) . (13) 

It is already obvious from the lower bound 1.2 on p 2 we have found for such wave functions 
eq . fliTf ) . In fact, choosing /3 variationally, we estimate (3 ~ 0.57 for G10- From this (3, one 
would find p 2 ~ 1.2 which is sensibly higher than the exact value (+ 17 %). 

The stability of the result p 2 ~ 1 can be understood from the following interest- 
ing fact concerning the relativistic wave functions. The integral giving p 2 from the 
momentum space wave function comes essentially from the intervall : 0.4 GeV < p < 2.4 
GeV . We find that in this interval, for the three relativistic models considered, the 
wave functions are rather close to each other (see fig. 1) and that they are rather well 
reproduced by an exponential form in r : 

if(r) = 2a~ 3 ^ 2 exp(— r/a) 



7 not the ISGW model for form factors 

8 Let us emphasize that we borrow from these authors only their spectroscopic model and not the 
calculation of form factors, which we do in our own way, a la Bakamjian-Thomas. 

9 In their paper, a s oc r, but anyway, the possible potential singularity would be smoothed out by the 
relativistic " smearing" 

10 Note that this value is quite different from the one quoted for GI in p5| Jlq] EM . 
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■m - (14) 

with a -1 ~ 0.75 GeV. Now it is found that for such values of a and m < 0.3 GeV, this 
exponential wave function gives p 2 ~ 1 (the contribution from outside the above interval 
is also small for this wave function). Two other remarks concerning the behavior in eq. 



(14) are in order. First the value a ~ 0.75 GeV does not correspond to the expected 
asymptotic behavior exp(— mr) ; we find that this latter behavior is only observed at 
very large distance r > 15 GeV -1 . Second, our value is quite in the range observed 



for wave functions of lattice NRQCD J2I2| when tuq — > oo : correcting the fit of [T^ 
for the lattice unit 1.75 GeV 0, we find that the lattice wave function is quite well 
reproduced for 0.5 < r < 5 GeV -1 by eq. (|14"D , with a -1 = 0.7 GeV. On the other hand, 
it must be emphasized that these relativistic wave functions are quite different from the 
non-relativistic ones, corresponding to the same mass spectrum. In particular the mean 



momentum y (0\p 2 |0) = 0.72 GeV for GI, is much larger than in the non relativistic ISGW 
model, where the same quantity is 0.5 GeV for the B meson ; the large momentum tail is 
also much higher. The various wave functions considered here are compared graphically 
in Fig.l. 




This correction is made with agreement of S. Veseli 
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Fig.l. Comparison of wave functions in momentum space multiplied by p ; 
the integral on p of the squared functions is 1 ; the unit on horizontal axis 
is 1 GeV ; solid line is for the GI spectroscopic model [17| ; long dashed line 
for the CNP model |2C 



; medium dashed line for the VD model [19| ; these 
three relativistic models are seen to give very similar wave functions ; short 
dashed line is for the non-relativistic ISGW model [14|, which is seen to yield 
a much shorter range wave function. 



Finally, let us note that the light quark mass m is not well determined in the above 
relativistic spectroscopic models. But at the same time p 2 is not very sensitive to m in 
the allowed range (e.g. m < 0.22 GeV for GI). This is due to the fact that m quite often 
enters through the combination y/p 2 + m 2 and that m <C \p\, which corresponds to the 
fact that the light quark is ultrarelativistic. 

Discussion of the effects of Relativity on p 2 

With the wave functions at hand, we can discuss more quantitatively the various 
relativistic effects affecting p 2 , which have been discussed in a qualitative manner in || 
using a formal v/c expansion in the light-quark velocity. We can first compare with the 
non-relativistic expression : 

p 2 = -m 2 {0\f 2 \0) . (15) 
3 

With such an expression, a non-relativistic spectroscopic model like ISGW would have 
given a very low value p 2 = .33 q completely unacceptable. A relativistic spectroscopic 



model of wave functions combined with the non-relativistic expression eq. flTq), if this 
makes sense, would have given a still lower and worse value, p 2 = 0.08 for GI, since the 
relativistic (0|f 2 |0) is smaller (e.g. for GI (0|f 2 |0) = 4.9 GeV 2 versus ISGW (0|f 2 |0) = 
9.1 GeV 2 ). 

Now, if we use, instead of eq. (|I5]), the expression of p 2 obtained from the relativistic 
treatment of the center-of-mass motion of B-T, and the relativistic free-quark current, 
eqs. (|5], |7|), the results are completely changed. These new relativistic effects (see the 
general discussion in ref. (||) shift p 2 upwards by a very large amount (around +1). This 
is particularly due to the fact that p 2 ce (second term in the r.h.s. of eq.(|7])) is increased 
by the momentum dependent factors po with respect to the non-relativistic expression 
(|j~5l). The shift (~ +1) of the relativistic expression of p 2 in terms of given rest frame 
wave functions, eqs. (§, 0), with respect to the non-relativistic expression, eq. (|1~5|), is 
in fact much greater than the value of the latter itself. This manifests the fact that the 



v/c internal light-quark velocity is not the least small (\J (0\p 2 |0) = .72 GeV for GI vs. 
m = 0.22 GeV) . For the non-relativistic spectroscopic model of ISGW, the result is 
now so large that it is again unacceptable. For a relativistic spectroscopic model like 
GI, the result is appreciably lower, so that it becomes reasonably compatible with QCD 
fundamental methods. 



12 This is the true value, given by the spectroscopic model itself. In the paper DM, radii are multiplied 
by a factor 2 to obtain a better fit to empirical radii. 
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Radiative corrections and comparison with QCD fundamental methods 



Although it could seem more natural to confront our results first with experimental 
data, we think that it is in fact a much more involved task than usually believed. In fact, 
comparison with other theories, although also requiring a careful discussion, is simpler 
and more conclusive, therefore we begin with it. One must warn the reader that many 
definitions of p 2 will be used in the following two sections, and that indeed some confusion 
has arisen about this question in the discussion of quark models. When comparing predic- 
tions of quark models with the ones of fundamental QCD methods like numerical lattice 
QCD or QCD sum rules, one always encounters a basic problem : the large QCD radiative 
corrections are obviously not included in the quark models. In particular quark models 
cannot describe the vhq — > oo limit of physical form factors, since these should have a 
logarithmic behavior (log rriQ) 1 due to the radiative corrections. These radiative correc- 
tions also manifest themselves in the HQET (heavy quark effective theory) perturbative 
evolution of the renormalized slope as function of the renormalisation scale /i. Using the 
relation given in|23|], and neglecting small a s terms [3 one has the general relation : 



81 a s {fi) 



which allows to define also a renormalisation-group invariant slope : 



Pren.gr.i™. = P 2 (P) + — ln "s(/i) • (17) 



16 

51 

Then, the standard assumption is that quark models describe the [i dependent 
matrix elements, obtained by factorization of the logarithmic coefficients Ci(mQ/p), at 
some particular low value of p, Pq. In this procedure, it is implicit that, at any other 
scale, we shall evolve p 2 (p) according to QCD (more exactly according to HQET in the 
present case) , which amounts to adding the QCD radiative corrections to the brute quark 
model result. p could not be much larger than 1 GeV, since only light-quark degrees 
of freedom and soft momenta are involved in the vhq — > oo limit of the quark models 
for £. On the other hand, if we want perturbation theory to make sense, we cannot go 
much lower than a p such that a s (p) ~ 1. This sets bounds for p . With such bounds, 
we can establish more precisely how to compare the renormalization-group invariant 
Pren.gr.inv. °f QCD fundamental methods with the p 2 of quark models. What we must do 
is to evaluate in the QCD fundamental methods p 2 (po) which is directly comparable with 
quark models. 

We start from the range of values of the renormalization-group invariant slope 
Pren.gr.inv. > predicted by fundamental QCD methods, which we denote by an additional 
QCD subindex. Lattice QCD has very large errors, which by themselves would give a very 
large allowed range 0.4 < Pq CD ^ ren .gr.inv. ~ 1-5, compatible with almost any estimate. 
However, taking into account QCD sum rules, we retain the narrower range : 

0-5 ^ PQCD, ren.gr.inv. ~ 1 ■ (18) 

Using eq.flTTI) : 

16 

Plen.gr.inv. = PVo) + gj hltt s (/i ) , (19) 



13 Here and in the following, only leading logarithmic corrections will be given for sake of simplicity, 
since we aim only at illustrating the effect of radiative corrections 
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with the range of po indicated above, this leads to : 

0-5 < pIcdM < 1.2 • (20) 

One can conclude that, as announced in the introduction, our quark model pre- 
diction /) 2 k 1, which is to be compared with Pqcd(Po), is quite compatible with the 
QCD methods, although closer to the upper edge of the range. A value like 1.4 |7|] seems 
improbable. 



Comparison with experiment 



Comparison with experiment requires a rather detailed discussion. There are two 
aspects which complicate the comparison : the first one is in fact mainly theoretical and 



has been already extensively discussed 2A\ the second one is purely experimental and has 
been perhaps underestimated. 

1) The first aspect is that the slopes at g^ aa . of the physical quantities that are 
measured in B —> D*tv are not identical with the p to be predicted by quark models. 
There are three reasons for that : 

a) the scale corresponding to the physical form factors may be taken around m c , 
then, considering the slope of the physical form factor , called p\ , the difference with 
the preceding scale p , corresponding to the quark model predictions, implies, grossly 
speaking : 

p 2 K) = P 2 (^o) - ^ In^T • (21) 
81 a s (p ) 

This evolution is of course due once more to the radiative corrections, and is based on 
the assumption that we must add these corrections to the quark model brute prediction, 
which of course cannot include them. 

(3) the physical form factor is also affected by I/tuq corrections, which are not 
included in our model, formulated at rriQ = oo. From the discussion of Neubert [24], one 
concludes that there is a large theoretical uncertainty on these corrections. Nevertheless, 
from the estimate of the U N^ 1 V ratio, Table 5.1. of Neubert, p. 372 , Na 1 ~ 1, one may 
expect that the 1/m corrections largely compensate numerically the radiative corrections 
in such a way that finally : 

2 2 
P A\ ^ Pren.gr.inv 

(=p 2 (//o) + i lna fl (/x )) • (22) 

7) some experiments do not measure directly p\ which is the logarithmic slope of 
h^, the main form factor involved in B — > D*£i>, but a quantity p 2 related to the slope 
of dT/dw. According to Neubert, p 2 can be safely related to p\ according to : 

p 2 «/4- 0.22 . (23) 

where 0.22 represents 0(1 /tuq) corrections due to the contribution of form factors other 
than h^ 1 - Thro ugh the above equations, one can a priori deduce straightforwardly p 2 (po) 
from the experimental measurements of p 2 or p\ , and compare with the quark model for 
the range of /i defined above. 
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2) However, there is still a second aspect, perhaps the main one, which makes the 
discussion more complex : the experiments, in view of the small number of experimental 
points, do not actually measure slopes at w m i n = 1, but rather average slopes over 



the whole range of w, 1 < w < 1.5 [26]. One can make fits which determine the slope 
at the origin if one fixes the functional form of the w dependence by a one parameter 
function (linear, exponential, ...), but the value of the slope at the origin very strongly 
depends on the chosen function : if the curvature is positive, the true slope may be 
much higher than the usually quoted numbers, which correspond to a linear fit, 
i.e. to an average slope ; moreover, the functional form is not strongly constrained by 
the data ; therefore, the true zero-recoil slope is not well determined. This phenomenon 
is generally underestimated when one compares experiment and theory. In our opinion, 
the simplest thing to be done is then to compare theoretical and experimental average 
slopes (average slopes will be denoted by bars). Theoretically, we define such an average 
slope by calculating the slope of the straight line 

iun(w) = 1 - p 2 (w - 1) (24) 

such that : 

Ei£K)-&«K)l 2 (25) 

i 

is minimum for a sufficient large set of points i in the interval 1 < Wi < 1.5. 

The result obtained in this way when one uses the GI model for calculating the 
wave functions is : 

p 2 GI = 0.75 . (26) 

It is sizeably lower than the number eq. (0), because the curvature of £ in our calculation 
is positive 0. In fact, our £ function is approximately given for the GI case by: 

/ 2 \ 2 ^ 2 

« M (i+^) • < 27 > 



As to the experimental numbers, we consider the value of p\ [26], since it is the 
one which is the most directly related to the theoretical calculation. The value obtained 
in the linear fit is in fact to be identified with p\ , according to eq. (|24]), and we quote 
therefore : 

p\ = 0.91 ±0.15 ±0.06 . (28) 

As we have explained, the quark model p 2 is to be identified with p 2 (p, ). Taking 
into account the relation between p 2 (p>o) and p 2 en , gr , inv ,, the relation between ^ ren .gr.inv. 
and Ka x (Iia 1 ~ Cren.gr.inv.) (see eq. |2~2|) , one should expect, according to the value of p,o : 

pl~P 2 Mtop 2 M + 0.2 . (29) 



Therefore, the comparison between eq. ( P6|) and eq. fl28|) is quite encouraging. 

However, the corresponding result for the ISGW spectroscopic model case is 
pjscw = 0-95 as calculated by us ; this is still compatible with experiment. Therefore, 
experiment would certainly not exclude presently the latter spectroscopic modelQ. It 

14 Note that our theoretical results for curvature are not far from the relations proposed by |p5| 
15 Let us once more insist that we are not speaking of the ISGW model of form factors, but of their 
spectroscopic model, when inserted in the B-T formalism for form factors 
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seems excluded only by comparison with the QCD methods. A conclusion is that, due 
to the phenomenon of the curvature of £, the present experiments are not testing 
accurately the slope at the origin, and they are not discriminating strongly between 
models through the average slope p 2 . 

Finally, let us note that a similar phenomenon occur at a purely theoretical level, 
in the analysis of UKQCD lattice results. The slope p 2 is estimated through a fit, and 
the result again depends on the functional form of the fit. Happily, the dependence is 
much weaker. The authors have chosen a particular form to get the result eq. @. A 
quadratic fit with unconstrained curvature gives instead 1 as central value ; according to 
us, this latter value is more significant as long as one has no firm certitude as regards the 
magnitude of curvature. 

In addition, as regards comparison of fundamental QCD methods with experiment, 
it is worth mentioning that in that case also, one should compare experiment with the 
average slope rather than with the higher slope at origin. 

Conclusion 

In conclusion, we have shown that the relativistic quark models of form factors a la 
Bakamjian-Thomas, combined with a relativistic spectroscopic model (i.e. with a square 
root kinetic energy) to calculate the needed wave functions, yield p 2 m 1. The chosen 
spectroscopic model was primarily the Godfrey-Isgur model, but the conclusion was also 
found valid for other similar models. This result contrasts with results obtained from 
a non-relativistic formulation of form factors and (or) a non-relativistic spectroscopic 
models: in the various cases, these results are either much too low or much too high. 
Our result p 2 ~ 1 is manifestly in better agreement with QCD fundamental methods, as 
illustrated by Table 1. 
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Quark models 


Form factors : 


Form factors : 




NR treatment 


relativistic treatment (B-T) 


NR spectroscopic 
model : ISGW |14j 


p 2 = 0.33 


p 2 = 1.283 






p 2 = 1.4 






p 2 = 1.64 § 


Relativistic 






spectroscopic models 






GIQ 


p 2 = 0.08 


p 2 = 1.026 ** 


GI 




p 2 = 1.15 ** 


vd 


p 2 = 0.124 


p 2 = o.98 ** 


cnp 


p 2 = 0.026 


p2 = 0.97 ** 


QCD fundamental 






methods(^ enflr i n „) 






QCD sum rules 


p 2 = 0.55 to 1 


Lattice QCD p 


p 2 = 0.9 +2 ( s tat) + ^ (syst) 
p2 = 1.2+2( stat) +2 (syst) 



Table 1. Theoretical values of p 2 : quark models, QCD sum rules, lattice QCD. 
NR stands for non-relativistic. The starred results are the quark model results 
advocated by us, with both relativistic wave equation and relativistic form factor 
formulation. The table illustrates the improvement and stability gained when including 

the two features together. All the results are explained in detail in the text. The 
numbers for NR form factor treatment are taken for conveniency at the B, but should be 
very close to the true infinite mass limit. The numbers for the relativistic treatment are 
taken at very large b mass. For quark models, citations refer only to spectroscopic quark 
models and the p 2 are the results of our own calculations, except when the reference is 

quoted after the value of p 2 . 
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